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Abstract. In cores of supernovae and crusts of neutron stars, nuclei can adopt interesting shapes, 
such as rods or slabs, etc., which are referred to as nuclear "pasta." In recent years, we have studied 
the pasta phases focusing on their dynamical aspects with a quantum molecular dynamic (QMD) 
approach. In this article, we review these works. We also focus on the treatment of the Coulomb 
interaction. 



I. INTRODUCTION 

In ordinary matter on the earth, atomic nuclei are roughly spherical. This fact can be 
understood in the liquid drop picture of the nucleus as being a result of the surface 
tension of nuclear matter, which favors a spherical nucleus, being greater than the 
Coulomb repulsion between protons, which tends to make the nucleus deform. However, 
in supernova cores and neutron stars, the situation changes completely (see, e.g., Refs. 
IH 12]): when the density of matter approaches that of atomic nuclei, i.e., the normal 
nuclear density po, nuclei are closely packed and the effect of the electrostatic energy 
becomes comparable to that of the surface energy. Consequently, around a density 
p ~ po/2, the energetically favorable configuration could be rod-like or slab-like nuclei 
embedded in the gas phase, or rod-like or spherical bubbles of the gas phase embedded in 
nuclear matter 0,|4|]. Such phases with exotic shapes of nuclei are referred to as nuclear 
"pasta" phases. 

Properties of the E^sta Erases in equilibrium states have been investigated by many 
authors so far 0.13. Isl BUS- These earlier works have confirmed that, for various 
nuclear models, the nuclear shape changes in the sequence sphere, cylinder, slab, cylin- 
drical hole, spherical hole, uniform with increasing density. This conclusion holds for 
the case of non-zero temperatures with a constant entropy Jql as well as for the zero 
temperature case. However, in these works (except for Ref. [5]) they took account of 
only several specific nuclear shapes and determined the energetically favorable one by 
calculating the energy for these assumed structures within a liquid drop model or the 
Thomas-Fermi approximation. Thus the phase diagram at subnuclear densities and the 
existence of the pasta phases should be examined without assuming the nuclear shape. It 
is also noted that at typical temperatures of the collapsing cores, several MeV, effects of 
thermal fluctuations on the nucleon distribution are significant. However, these thermal 
fluctuations cannot be described properly by mean-field theories such as the Thomas- 



Fermi and the Hartree-Fock approximations. 

In the processes of supernova explosion and succeeding neutron star formation, the 
pasta phases can be formed in two stages. During the gravitational collapse, matter in 
the collapsing core is compressed and the central region starts to solidify into a bcc 
lattice because the Coulomb repulsion among nuclei gets stronger. In the later stage 
of the collapse, the central density reaches po and the pasta phases could be formed 
from the bcc lattice due to the compression. After a bounce of the core takes place, the 
temperature of the core increases to O(10) MeV and nuclei melt completely. As the 
cooling process of the protoneutron star proceeds later on, the pasta phases as well as 
normal spherical nuclei could be produced again from hot uniform nuclear matter at 
subnuclear densities. 

The previous studies, which are based on a comparison of the energy with an as- 
sumption of the nuclear shape, cannnot answer whether or not the pasta phases can be 
formed dynamically within the time scale of the neutron star cooling nor whether or 
not the transitions between them, which is accompanied by drastic changes of nuclear 
shape, can be realized under nonequilibrium conditions in the collapsing cores. To solve 
these problems, molecular dynamic methods for nucleon many-body systems are suit- 
able. They treat the motion of the nucleonic degrees of freedom and can describe thermal 
fluctuations and many-body correlations beyond the mean-field level. 

Using one of the molecular dynamic methods, the quantum molecular d yna mics 
(QMD) full , we have arrived at the two following major conclusions lfl2l[l3lTl4tl . (1) 
The pasta phases are formed from hot uniform nuclear matter by cooling it down within 
a time scale of ~ 0(1O 3 — 10 4 ) fm/c. This supports the idea that the pasta phases exist 
in neutron star crusts. (2) The transition from rod-like nuclei to slab-like nuclei and that 
from slab-like nuclei to rod-like bubbles can be realized by compression of matter. This 
suggests that the pasta phases could be formed in collapsing cores. 

Throughout the present article, we set the Boltzmann constant &b = 1 ■ 



II. METHOD: QUANTUM MOLECULAR DYNAMICS 



Among various versions of the molecular dynamic models, QMD 111 111 is the most 
practical and suitable for studying the pasta phases. This method is so efficient that 
we can describe rod-like and slab-like nuclei in terms of nucleon degrees of freedom 
even though they consist of macroscopic number of nucleons. In addition, we focus on 
the pasta structure and the configuration of nuclei on a macroscopic scale, where the 
shell effects lll5l ll6l Intl . which cannot be described by QMD, may be less important 
(in the case of supernova matter, shell effects are always small). Thus QMD is a good 
approximation for our present problem. 



1. Model Hamiltonian 



In our studies, we use a nuclear force given by a QMD model Hamiltonian with 
the medium-equation-of- state parameter set in Ref. 1X811 . This Hamiltonian contains 



the momentum-dependent "Pauli potential" to reproduce the effects of the Pauli prin- 
ciple phenomenologically. The Pauli potential generates repulsive forces between two 
identical particles close together in phase space. Parameters in the Pauli potential are 
determined by fitting the kinetic energy of the free Fermi gas at zero temperature. 

This Hamiltonian reproduces the binding energy of symmetric nuclear matter, 16 
MeV per nucleon, at the normal nuclear density po = 0.165 fm~ 3 and other saturation 
properties: the incompressibility is set to be 280 MeV and the symmetry energy is 34.6 
MevOl. This Hamiltonian also well reproduces the properties of stable nuclei relevant 
for the situations we consider: the binding energy (except for light nuclei from 12 C 
to 20 Ne) JH, and the rms radius of the ground state of heavy nuclei with A > 100 
lfl9ll . It is also confirmed that another QMD Hamiltonian close to this model provides a 
good description of nuclear reactions including the low energy region (several MeV per 
nucleon) [20], which is one of the essential elements when one studies the dynamical 
processes as in Sect. III-2. 



2. Coulomb Interaction and its Screening Effect 

The Coulomb interaction is one of the essential elements for the pasta phases. One 
must carefully treat the long-range nature of the Coulomb interaction since electron 
screening is negligibly small in this system fllzHIzS]- A key quantity in determining 
the importance of electron screening is the ratio of the scale of the inhomogeneity, 
typically half the internuclear spacing R, to the Thomas-Fermi screening length Atf = 
[Ane 2 {dnf^ /djlf^) (o)]~ 1//2 = yJn/Anikf"') -1 , where ny and jif^ are the averaged 



electron number density and chemical potential, respectively, and — y 3n 2 n^ 
is the electron Fermi wave vector. We note that, in the density region of the pasta 
phases, R/Xj-p < 1 lE3l flll. which means that it is a good approximation to assume that 
the electron density distribution is uniform and neglect the screening effect; but one 
must take account of the long-range nature of the Coulomb force by, e.g., the Ewald 
summation procedure II 1211 . 

It is fortunate that the screening effect is small; otherwise one has to solve equa- 
tions of motion for the electronic degrees of freedom together with that for the nucle- 
onic ones because the coupling between these two components cannot be neglected. 
Thus introducing a screening length shorter than the actual value for calculating the 
Coulomb force between protons keeping the density of electrons uniform is incon- 
sistent. Let us restrict ourselves to the equilibrium state and discuss this point more 
quantitatively using a liquid drop model in the Wigner-Seitz approximation. The energy 
density of a relativistic degenerate electron gas with an inhomogeneous density distri- 
bution n e (r) = + 8n e (r) is given by E e [n e (r)] = E^\n^) + 8E e ~ Ee°\ne)[l + 

(2/9V c ) .Lu d 3 r (5n e (r)/4 0) ) 2 ], with E^\n { e ] ) = (3 / 'A)hckf' ] nf ] ', where V c is the vol- 
ume of the Wigner-Seitz cell. If the screening effect is significant, the Coulomb energy 

among protons decreases, but on the other hand, the (8n e (r) /n^) 2 term cannot be ne- 
glected. Calculations assuming a small Atf and uniform n e artificially underestimate the 
Coulomb energy among protons but discard the significant contribution of 8E e . Now we 



estimate 8E e for an unrealistically small value of Atf- We take the phase with spherical 
nuclei in the density region where the pasta phases start to appear, i.e., the volume frac- 
tion of nuclei is u — (r^/r c ) 3 ~ 1/8 according to the condition of the fission instability 
of spherical nuclei [1]. Here and r c are the radii of nuclei and the Wigner-Seitz cell, 
respectively. Taking supernova matter for simplicity, in which the density of dripped 
neutrons is negligible, with the nucleon density in nuclei n = 0.1 fm 3 , proton fraction 
in nuclei *n = 0.3, r c = 15 fm, and = 15/2 fm, and using expressions in Ref. lElll . 
8E e per nucleon, 8E p /uk, is estimated as 8E e /n^ ~ 0.23 MeV 1 for Atf = 10 fm (this 
value is adopted in Ref. [24]). Note that this is much larger than the energy difference 
between different pasta phases in the same density region in supernova matter, O(10) 
keV per nucleon. 



III. SIMULATIONS AND RESULTS 

In the present section, let us review our previous works iflU 13, [3]. In our QMD 
simulations, we treated a system which consists of neutrons, protons, and electrons in a 
cubic box with periodic boundary conditions. The system is not magnetically polarized, 
i.e., it contains equal numbers of protons (and neutrons) with spin up and spin down. 
The relativistic degenerate electrons which ensure charge neutrality are regarded as a 
uniform background The Coulomb interaction is calculated by the Ewald 

method taking account of the Gaussian charge distribution of the proton wave packets. 



1. Realization of the Pasta Phases and Equilibrium Phase Diagrams 

In Refs. O0, 

we showed that the pasta phases are produced from hot uniform 
nuclear matter by cooling it down and we studied phase diagrams at zero and non- 
zero temperatures. In these works, we first prepared a uniform hot nucleon gas at a 
temperature T ~ 20 MeV. We then cooled it down slowly until the temperature got ~ 0. 1 
MeV or less for O(10 3 — 10 4 ) fm/c, keeping the nucleon number density constant. Note 
that this cooling time scale is much larger than the time scale T re i ax for relaxation of 
the system, which is estimated as T re i a x ~ (length scale of the inhomogeneity)/(sound 
velocity) ~ 10 fm/0.1 c = 100 fm/c. In the cooling process, we mainly used the 
factional relaxation method, which is given by the QMD equations of motion plus small 
friction terms: R, = dJ^/dPi - % R dJf/dRi, P ; = -dJf/dRi - ^pdJ^/dPi, where 
is the QMD Hamiltonian, R ; and P ; are the position and momentum of nucleon i, 
respectively, and > are the frictional coefficients. 

The resulting nucleon distributions of cold matter at x = 0.5 are shown in Fig. [T] 
We see from these figures that the phases with rod-like and slab-like nuclei, cylindrical 
and spherical bubbles, in addition to the phase with spherical nuclei are reproduced. 



This calculation is based on the linearized Thomas-Fermi approximation and thus the true value of 
8E e /riN would be larger than this value. 



FIGURE 1. Nucleon distributions of the pasta phases in cold matter at x — 0.5; (a) sphere phase, 
0.1 po; (b) cylinder phase, 0.225po; (c) slab phase, 0.4po; (d) cylindrical hole phase, 0.5po; (e) spherical 
hole phase, 0.6po- The red (darker) particles represent protons and the green (brighter) ones neutrons. 
Taken from Ref . II 211 . 



The above simulations show that the pasta phases can be formed dynamically from hot 
uniform matter within a time scale ~ O(10 3 — 10 4 ) fm/c. 

In Fig. El we show snapshots of the nucleon distributions for p = 0.225po at T = 1,2 
and 3 MeV. This density corresponds to the phase with rod-like nuclei at T = 0. We 
have observed the following qualitative features: at T ~ 1.5 — 2 MeV the number of 
evaporated nucleons becomes significant; at T > 3 MeV, nuclei almost melt and the 
nuclear surface is hard to identify. 

The phase diagram for x = 0.5 is plotted in Fig. |3j The critical temperature of this 
model is > 6 MeV. In the region below the dotted lines at T < 3 MeV, where we can 
identify the nuclear surface, we have obtained the pasta phases in the same sequence as 
in the earlier works: from lower densities, spherical nuclei [region (a)], rod-like nuclei 
[region (b)], slab-like nuclei [region (d)], cylindrical holes [region (f)], and spherical 
holes [region (g)]. In addition to these pasta phases, structures with multiply connected 
nuclear and bubble regions (i.e., sponge-like structure) such as branching rod-like nuclei, 
perforated slabs and branching bubbles, etc., have been obtained in the regions (c) and 
(e). Further study using a larger system is necessary to conclude the existence of these 
sponge-like phases (see also Refs. II lL l25L 12610 . 



2. Structural Transitions between the Pasta Phases 

In Ref. JH, we demonstrated the second point mentioned at the beginning of this ar- 
ticle. We performed QMD simulations of the compression of dense matter and observed 



(a) T=1 MeV (b) T=2MeV (c) T=3MeV 




FIGURE 2. Nucleon distributions for x = 0.5, p = 0.225p at temperatures of 1 , 2 and 3 MeV. The total 
number of nucleons N = 16384 and the box size Lbox =76.14 fm. The upper panels show top views along 
the axis of the cylindrical nuclei at T — 0, and the lower ones show side views. Protons are represented by 
the red (darker) particles, and neutrons by the green (lighter) ones. Taken from Ref. 
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FIGURE 3. Phase diagram of matter at x = 0.5 plotted in the p - T plane. The dashed and the dotted 
lines on the diagram show the phase separation line and the limit below which the nuclear surface 
can be identified, respectively. The dash-dotted lines are the phase boundaries between the different 
nuclear shapes. The symbols SP, C, S, CH, SH, U stand for nuclear shapes, i.e., sphere, cylinder, slab, 
cylindrical hole, spherical hole and uniform, respectively. The parentheses (A,B) denote an intermediate 
phase between A and B-phases with a multiply connected structure. Simulations have been carried out at 
points denoted by circles. Adapted from Ref. 1 13]. 
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FIGURE 4. Snapshots of the transition process from the phase with rod-like nuclei to the phase with 
slab-like nuclei. The red (darker) particles show protons and the green (lighter) ones neutrons. The box 
size is rescaled to be equal in this figure. Adapted from Ref. (3- 



the transitions from rod-like nuclei to slab-like ones and from slab-like nuclei to rod-like 
bubbles. 

The initial conditions of the simulations are samples of the phase with rod-like nuclei 
(p = 0.225po) and of the phase with slab-like nuclei (p = 0.4po) with 16384 nucleons 
at x = 0.5 and T ~ 1 MeV. We adiabatically compressed these samples by increasing 
the density to the value corresponding to the next pasta phase taking O(10 4 ) fm/c. 
This time scale is much larger than the typical time scale of the deformation and 
structural transition of nuclei (e.g., that of nuclear fission is ~ 1000 fm/c). Therefore, 
the compression in the simulations is adiabatic with respect to the change of the nuclear 
structure, so that the dynamics of the structural transition of nuclei observed in the 
simulations is physically meaningful, and is essentially independent of the compression 
rate, etc. 

The resulting time evolution of the nucleon distribution is shown in Figs. |4| and In 
Fig.|4j we see that the phase with slab-like nuclei is finally formed [Fig. |4]-(6)] from the 
phase with rod-like nuclei [Fig. |4]-(1 )] . We note that the transition is triggered by thermal 
fluctuations, not by the fission instability: when the internuclear spacing becomes small 
enough and once some pair of neighboring rod-like nuclei touch due to the thermal 
fluctuations, they fuse [Figs. |4[(2) and|4[(3)]. Then such connected pairs of rod-like 
nuclei further touch and fuse with neighboring nuclei in the same lattice plane like a 
chain reaction [Fig. El- (4)]; the time scale of the each fusion process is O(10 2 ) fm/c. 

The transition from the phase with slab-like nuclei to the phase with cylindrical holes 
is shown in Fig. |5j When the internuclear spacing decreases enough, neighboring slab- 
like nuclei touch due to the thermal fluctuation as in the above case. Once nuclei begin to 




FIGURE 5. The same as Fig.|4]fbr the transition from the phase with slab-like nuclei to the phase with 
cylindrical holes (the box size is not rescaled in this figure). Adapted from Ref. Q. 



touch [Fig. El-(2)] , bridges between the slabs are formed at many places on a time scale 
of O(\0 2 ) fm/c. Initially, the bridges cross the slabs almost orthogonally [Fig. I3-(3)]. 
Nucleons in the slabs continuously flow into the bridges and the bridges become wider 
and merge together to form cylindrical holes. Then the cylindrical holes finally relax 
into a triangular lattice [Fig.|5[(6)]. 



IV. SUMMARY AND FUTURE PROSPECTS 

Using QMD, we have shown the following two things: formation of the pasta phases in 
neutron star crusts by cooling and formation of rod-like bubbles from slab-like nuclei 
and that of slab-like nuclei from rod-like ones in supernova cores by compression of 
matter. In closing, we list important issues to be clarified in the future. 

1. Formation of the pasta phases by compression [27] 

A remaining problem is a transition from a bcc lattice of spherical nuclei to a trian- 
gular lattice of rod-like nuclei induced by compression. If this process is confirmed, 
existence of the pasta phases in supernova cores will be almost established. 

2. Effects of uncertainties of nuclear forces & properties of nuclei 

We employed a specific nuclear force. Although we have confirmed that this nu- 
clear interaction yields reasonable values for the surface energy and the proton and 
neutron chemical potentials in neutron matter, etc. Ill 211 . in addition to reproducing 
saturation properties [18], further systematic survey is needed to understand effects 
of uncertainties in nuclear properties on the pasta phases Il28ll . Especially, uncer- 



tainties in the surface properties are the most important elements to be examined. 
3. Detailed study of astrophysical consequences 

If formation of the pasta phases in actual astrophysical situations is established, 
quantitative discussion about effects of the pasta phases on astrophysical phenom- 



ena will be increasingly important. As we have stressed previously lllflLll2ll . effects 
of the pasta structure on coherent scatterin g of neutrinos in the stage of the neutrino 
trapping is an interesting problem l!24ll3Ll28Ll29ll . 
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